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The   movement   of   sand   in  the   swash-zone   south  of    the 
head  of  Monterey   Canyon  was   studied   during  February  and 
March,    1966.      A  stationary   sampler   was   designed  and  used  in 
conjunction  with  dyed  fluorescent  sands   to  trace  the  rate  and 
direction  of  natural  sand  movement.      A  sequential   multiple 
linear  regression  program  was  used  to   statistically  analyze  the 
effects   of   this   canyon  and  several  other   environmental  para- 
meters on  the   movement  of  beach  sand.      In  all  observations 
made,    the  sand  was  found  to   move  toward  the  canyon  head. 
The   canyon  also  appears  to  be  a   major  factor  affecting  the 
rate  of  beach  sand  drift. 
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1.       INTRODUCTION 

The  Problem 

Monterey   Bay  provides   a   unique  area  for   the  investigation 
of   submarine  canyon  effects   on  beach-sand  migration.      The 
object   of  this   study  is   to   determine,    through  a   statistical 
approach,    the   effects   of    the  Monterey  Canyon,    as   well  as    that 
of   several   other   measured   parameters,    on  the   rate  and  direction 
of  swash-zone  sand   movement. 

In  recent   years,    a   considerable   degree   of   progress   has 
been  achieved  in  the  field  of  beach-sand  dynamics  and  transport 
processes.      Of  particular  interest  is   the  interaction  occurring 
between  the  beach  and  environmental  factors.      Also  of  interest 
is   the  relative  contribution  that   the  individual  parameters   make 
to   the  net   movement  of  beach   sand.      Due  to  the  large  number 
of  agents  interacting  in  this   situation,    the  problem  is   complex 
and  difficult  to   resolve  without  employing  a  statistical  approach. 
It  is   proposed   then,    that  by   measuring   the   environmental 
parameters  and  applying  a   statistical  analysis,    a   meaningful  set 
of  relationships   can  be  obtained. 

Background 
Chamberlain  £  1J    investigated   the  littoral   sand  budget  for 
a   series   of   California  beaches   associated  with     submarine   canyons. 
In   each   of    these  regions   there  is   a   littoral   cell   containing   stable 
beaches   where  a  balance   exists  between  sediment  input  and  loss. 
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Typically,    each   cell  is   bounded  on   one   side  by  rocky  headlands 
■where   the  influx   of   sand  is   restricted.      Within  the   cell,    the 
sand  is   supplied  to    the  littoral   zone  by   coastal   streams   and 
rivers,    and  is   transported  along  the   coast.       The  presence   of 
the  canyon  head  acts   as   a   terminal  point   for   the   cell   and   inter- 
cepts  the  littoral   sand   which  is   then  transported  down   the 
canyon.      Previously,    Sheppard    \_8J    noted   this    effect  in  a   series 
of   detailed   measurements    made  in  La  .Jolla   Canyon.      He  found 
that  the  head  filled  at   a  rate  approximately   equal   to   the  average 
net   transport  rate   of    sand   for   the   entire  region,    and   that  it 
subsequently  slumped  into   the   canyon   through    the   mechanism   of 
turbidity  currents. 

In  general,    the  net   result   of   a   submarine   canyon  with  its 
head   extending  into   shallow   water   is  a   central   region   of   low 
wave   energy   surrounded  by   regions    of  higher   wave   energy  L1QJ  • 
In  Monterey  Bay,    refraction  patterns   clearly  indicate  a   region 
of   energy   divergence   located   directly  over   the   canyon  with 
adjacent    zones   of   energy   convergence      (Fig.    1).       Conservation 
of   energy   and   continuity   considerations   provide  for   a  net   move- 
ment   of   water   into    the  divergent   area,    resulting  in  average 
littoral   currents    directed   toward   the   canyon  from   both   sides. 
The  largest  instantaneous   rates   should   occur   nearest   the   canyon 
head    (_9j  •      This   line  of   reasoning   can  also  be   extended   to   include 
the  littoral   drift  problem,    where  again,    the  largest   rates   would 
be   expected  in  the  vicinity   of    the   canyon  head. 
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Some  indication  of   these   conditions   occurring  in  Monterey- 
Bay  has   been  demonstrated  in   two   separately   conducted   observa- 
tions.     A  study   made  by   Hohenstein,    Jaeger,    and  Jones    \jf\    on 
Del  Monte  Beach  at  the   southern   end   of    the  bay  provided   evidence 
of   little   or  no   sand   drift  occurring  on   the  beaches   in   that   region. 
A  trial  investigation   conducted  by  the  authors   in  January   1966, 
one   quarter   mile   south   of    the   canyon  head    (point  B    in  Fig.    2), 
indicated  a   strong  rate   of   drift   toward   the   canyon  head. 

These  previous   remarks   all   suggest   that   the   presence   of 
the  Monterey   Canyon   might  indeed  have   some   effect   on  beach- 
sand   movement.      Since  no   known  studies   relating  both   canyon 
and   environmental   effects    to   beach   sand   migration  have  been 
conducted,    a   detailed  investigation  of    this   problem    was    desirable. 

Organization 
The  first   part   of    this   thesis   is   concerned  with  the   approach 
to   the  problem,    including  techniques   available,    methods    employed, 
and  related  assumptions.      In   the   next  part,    the  procedures 
involved  in  data   collection,    data  analysis,    and   data  processing 
are   discussed.       Part   three   is   a   presentation  and  analysis   of   the 
processed   data,    while   part   four   is   a   discussion  and  interpretation 
of   these   results,    with  related   conclusions.      The  final   section 
deals    with   recommendations   for  future   related   studies. 
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FIGURE   2   -   STATION    GRID    LOCATIONS 
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2.       STUDY    TECHNIQUES 

Discussion 
The  basic   problem   of   sand   movement  is   a   three-dimensional 
study  in  time,    having   offshore,    longshore,    and  vertical   components 
A   simple  problem   involving  a   one-dimensional   rate   can  be  treated 
by  ignoring   the   offshore  and  vertical   components.      The   resulting 
problem    can  then  be  divided  into   two   easily  investigated   sub- 
problems.      The  first   is   a   determination   of    the   direction  of    sand 
movement   with   respect   to   the   canyon  head,    and   the   second  is 
a   determination  of    the  rate   of   sand   movement   as    related   to   the 
canyon   effects   in   conjunction   with   other   environmental   parameters. 

Sand  Movement 

In  the  past,    determination   of    sand   movement  has   been  a 
difficult  problem.      Various    techniques   have  included  the  use   of 
physical   characteristics    such  as    color,    shape,    size   distribution, 
mineral   content,    and   magnetic   properties   as   identifiable  features 
to   trace  sand   movement.      More  recent    methods   primarily   employ 
the  use   of   radioactive  particles   for   tracing  silts   and   clays, 
painted   or   dyed   tracers   for   pebbles   and   cobbles,    and  dyed  fluor- 
escent  grains   for   sands   £3J  • 

A  representative   tracer  particle  used  in   a   sand   movement 
study   must  react   to   fluid  forces   in  a   manner    closely  resembling 
that   of    the  naturally   occurring  grains,    while    maintaining  its 
physical   identity.      This   requires    the   tracer   to  have  a   size, 
distribution,    density,    shape,    surface   chemistry,    and   strength 
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similar   to   the   natural   beach   material.      In  theory  then,    if    most 
of   the  above   criteria  are   satisfied,    a   tracer   may  be  used  to 
determine  a  representative  value  for  both   the   direction  and   rate 
of   movement   of   the  naturally   occurring   sands.      Ideally,    the 
tracers   should  be  prepared  from   sands   in  the  beach  area   where 
the  study  is   to   be   conducted.      It    may  be  necessary,    however, 
to   prepare  artificial   tracers    or   natural   tracers   obtained  from   a 
different   locale.      Care   must  be   exercised  to    match   the  tracer 
characteristics    to   those  of   the   studied  area    [llj   • 

Sampling  and  analysis   techniques    may   vary   depending   on   the 
results   desired.      A  space   integration   sampling   method,    where 
tracer   counts   are  made   over  a   multi-sample  grid   system,    provides 
information  on   the  direction  and   distribution  of   sand   movements. 
Subsequent   sampling   without   replacement   of   the   tracer   can  be 
used   to   determine   a   time   rate   of   particle   distribution.      This 
method  is    slow   and  requires   analysis    of   a   large  number   of   sand 
samples.      A  faster,    less   tedious   procedure  is   the  time  integration 
method,    where  a  few    samples   are   obtained  in  a  line  normal  to 
the  littoral   drift.       Periodic    sampling  at   a  given   distance  from 
the  injection  point   will   establish  a   rate  by  providing   the   time  at 
which  a   significant   number   of    marked   grains   reaches   the   sampling 
point    [is]    .      Laboratory   determinations   in   conjunction   with   the 
above   sampling   methods    can  provide  both   qualitative  and   semi- 
quantitative results .       The   direction  and  rate  values   obtained  are 
qualitative  in  that   only  relative   concentrations    of    the   tracers 
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are  important.      Quantitative   observations   involve   major  assump- 
tions  as   to   the  vertical   and   seaward  distribution  of  the   tracer 
particles. 

Interactions 

Reliance  upon  information   supplied  by  direct   observations, 
in  which  inferences   are   made   regarding   the  physical   relationships 
of   the  agents   acting,    suggests    that   the   study   of   beach   processes 
should  be  primarily   statistical.      This   approach   provides   a   ration- 
al basis   for   the  study   of   beach  phenomena  because  it  includes 
the  following  features.      A  specific   set    of  variables   is    measured, 
and  a  plan  for   sampling   at   random   times   and   places   is   designated. 
An  analysis   of   the  data  locates   the   sources   of   variability  and 
the   results   are  used   to   predict   relationships   among  these 
agents    [6 J  . 

There  are  several  statistical   models   which   can  be  used. 
Multiple  regression  analysis,    correlation  analysis,    factor   analysis, 
and  discriminant   function  analysis   are  all  useful  in  particular 
situations   [_7]  • 
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3.       THE    STUDY 

Location 
The  field   phase   of   this   study  -was   conducted   on  a   stretch 
of  beach   extending   southward  from   the  Moss   Landing 
Pier    (Fig.    2).       Physical   sampling   complexities   restricted      the 
study   to   considering  only   those   movements   in   the   surface  layer 
of   the   swash-zone.       Time  and   transportation   deficiencies   also 
limited   the   study     to      one        side   of    the   canyon  and  a   longshore 
extent   of   approximately   one    mile.       The  head   of   the   canyon   was 
chosen   to   be  the  point   where   the   offshore   contours   sharply 
closed   the  beach   near   the   end   of   the   pier    (Fig.    2).      Inspection 
of   the   contours   one   mile   south   of    the   pier  indicated   that  this 
distance   would  be   sufficient   to   recognize   a   difference  in   canyon 
effects,    if   any   existed. 

Rate  and   Direction  of    Drift 
Since   the   study  area   was   a   sand  beach,    the   drift   rate  and 
direction  determinations   were   obtained  using   dyed  fluorescent 
tracers   in  a   time  integration   sampling   technique.       Preliminary 
studies    considered   only   sand   movements    toward   the   canyon  and 
sampling   was    conducted   on   the   canyon  side   of    the  injection  point 
(Fig.    3a).      Subsequent   evidence   revealed   that   single-side   sampling 
did  not  provide   discrimination  between   positive   drift  and   dispersion, 
Therefore,    a   single   sampler   with  two   opposing  injection  points, 
having   differently   colored   tracers,    was   used    (Fig.    3b).       Hence, 
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direction   of   drift  and  dispersion  could  be  positively  identified  by 
a   single  observation. 

Parameters    Used 


To    obtain  genuine   significance   from   this   study,    the   canyon 
effects    must  be   evaluated  in   conjunction  with   other    environmental 
factors   affecting   sand  drift.       In  all,    eleven  independent   variables 
■were  included  in  this   study. 

Distance  from   the   canyon  head   was    chosen  as   the  parameter 
to   represent    canyon   effects.      With  increasing   distance  from   the 
canyon,    decreasing   refraction   causes   a    more   equal   distribution   of 
wave   energy  on  the  beach,    thereby  reducing  the   divergent    effects 
of   the   canyon  and   the   expected  rate   of   drift.      Breaker   height, 
wave  period,    and  angle  of   approach   were   chosen   to   represent   the 
intial   energy  present   and  its    distribution  to   the   surf    zone. 
Height   of   the  tide   regulates   the  amount   of   the  initial   energy 
which   transits   the  surf    zone.       In   this   respect,    low    tide   causes 
a  large  amount   of    energy  to   be   expended  on  bottom   features 
prior   to   reaching   the   swash-zone.  |10J  .      Mean   slope  in  conjunction 
with   swash-length  and  period   was   used  as   a    measure  of   the 
amount   of    energy  being  released  to   the  beach.       The  latter   two 
were   combined  as   a   single  parameter,    swash-velocity,    to   provide 
an      indication   of   the  frequency  of   wave  action  on   the   swash-zone. 
The   mean   diameter   of   the  grains   that    moved   was   used   as 
an  indication  of   the   energy  being  transported  along   the  beach 
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since   large   grains   would   require   more   energy   than   smaller   grains 
to  be  transported  a   given  distance. 

The   -wind   was    considered  because   of   its   ability   to   locally 
alter   the   energy.       Two   factors   are  important:    the  longshore    wind 
component    can   contribute   to   the   lateral   distribution   of    energy; 
and  the  onshore/offshore   component   can   alter   existing   wave 
forms   fi+1  •      Longshore   current   with   its   ability  to    move   sand, 
has  been  classically  associated   with  the   energy  distribution  of 
the  aforementioned  parameters   and   was   therefore  included. 

Interaction  Analysis 
Sequential    multiple   linear   regression   was    the   statistical 
model   chosen   to  analyze  the   observed   interactions    |_7j  .       For 
this   approach   to   be   effective,    a   high   degree   of   variability   of 
the   observed  parameters    must   exist.      Ideally,    the  variable 
should  be   measured   over   the   maximum   range   of   values    which   can 
be   expected  during   the   study  period.      To   adequately  investigate 
the  area   and   still   provide   sufficient   variability  in  the   statistical 
approach,    a   minimum   of    eight   to   ten  grid  points,    with   sufficient 
grid   point   spacing,    is   required.      To    maximize   this    effect   an 
eight-point   grid   spaced  at   800-foot   intervals   was    established. 
To   provide   confidence  in  the   statistical   results,    a   minimum   of 
three   observations    per   variable   considered   was    required    {_4J   . 
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4.   PROCEDURES 

Data   Collection 

The  initial   step  in   the  data   collection  phase   of   the   study- 
was   the   establishment   of   a   sampling   station  grid.       This   was 
accomplished  by  using  a   transit   and   stadia   rod.      Stations    were 
laid   out  at   800   foot  intervals   using   the  Moss   Landing   Pier   as 
a  base  point    (Pig.    2). 

In   order   to   provide  an  accurate  indication   of   the  grain 
sizes   actually   moving   along   the  beach,    the   tracer   sands   used 
had  a   wide  grain  size   distribution.      Orange  and  green  fluorescent 
tracer   sands   ranging  in   size  from    0.065    to    1.0    millimeter   diameter 
were  utilized  for   this   study.       The   dyed   sands    were  an  artificial 
mixture   of   local   coastal   dune  and   sea   floor   sand   procured  by 
Hohenstein,    et  al,    for  a   1965   beach-sand  study   [5]  •      The 
mixture   was    coated  by   the   Great  American   Color   Company   of 
Los  Angeles,    California. 

During   the  period   of   2    February   through   22   March    1966, 
the   eight   sampling   stations    (Fig.    2)    were   each  visited  six   times 
on  a   random  basis.      Collections    were   made  under   as   widely 
varying   environmental   conditions   as   possible. 

The  first   seven  observations   used  fluorescent  blue  anth- 
racene tracers    seeded  in  the   swash-zone  at   the  grid   stations. 
At   a   position   50    feet   toward  the   canyon  from   the  grid   station, 
a   metal   can  cut   to   scoop   only   the  upper  half -inch    of   the  beach 
was   used  to    collect   a   sample  in  a  line   through   the   swash-zone. 
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A  different   method   of   using   two   colors   of   dyed   sand  and 
seeding   on  both   sides   of   a   stationary   sampler    (Fig.    2+)    was 
utilized   throughout   the   remaining  2+1    observation  periods. 

After      a       station  was    selected,    the   sampler   was    estab- 
lished in  the   mid   swash-zone.      It   was    secured  in  position  with 
four   two-foot   iron   spikes.       This   was   necessary  due  to   the 
occasional  violent  action   occurring  in  the   swash-zone. 

Prior   to   tracer   seeding,    a  five— minute   sample   was    obtained 
to   determine  a  background   count   of   dyed   sand   that   may  have 
been  present   from   previous   observations.       The   next   step   was 
to   seed  approximately  five   pounds    of   green  tracer    50    feet  from 
the   sampler  in  a   direction  away  from   the   canyon  while   simulta- 
neously seeding  an   equal  amount   of    orange   tracer   the   same 
distance   towards    the   canyon.       Seeding   was   accomplished  by 
injecting  the   marked   sand  in  a   line   normal   to   the  beach  immedi- 
ately behind  an  outgoing   swash    (Fig.    3)» 

To   preserve   sampling   continuity,    two    collection  boxes    were 
used  alternately  and   samples    taken  at   seeding   time   plus    5,    plus 
10,    plus   15,    plus  20,    plus  25,    plus  30,    plus  40,    plus   50,    and 
plus   60    minutes.      Samples    were   transferred  from   the   collection 
boxes    to   polyethylene  bags   and  appropriately  tagged  for   later 
analysis. 

During   the   course   of    each   sampling  period,    the   environmental 
parameters,    as   described  in  Appendix   I,    were   measured  and 
recorded. 
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Midway  through   the   data   collection  phase,    a   »Jeep   was 
obtained.       Use  of   this   vehicle  vastly   simplified  and  accelerated 
the  sampling  procedure  and   made  possible  a   greater   number   of 
observations  in  the  limited   time  frame  allocated  for   data 
collection. 

Data  Analysis 
After   field   collections,    the  36O    sand   samples   were  brought 
to   the  laboratory  for   analysis.       Each  sample   was   force-dried  in 
an  oven  at   250°  F.      It   should  be  noted   that   very  high  temp- 
eratures   ( 400-500°  F)    have  a   detrimental    effect   on  the  fluor- 
escent  properties   of   the   marked   sand,    hence   the  lower   drying 
temperatures   and   correspondingly  longer   drying   times   were  used. 
Following   the   drying   process   the  individual   samples   for    each 
observation  were  reduced   to   a   given  volume   to   provide  a    common 
reference.       The   samples    were  then   shaken  through  a   nest   of 
six   Tyler   sieves   having   mesh   openings   ranging  from    0.061   to 
1.981   millimeters.       Each  grain  size   was   inspected   individually 
under  an  ultra-violet   light  and  the  number   of    tracer   grains 
counted    (Appendix  II).      In   order   to   accelerate  and   simplify  the 
counting   of   tracer   grains   a   multiple-riffle  type   sample   splitter 

was   used.       The  number   of    marked  grains   in  a   sample   was   assumed 

s— 1 
to   be:    n=2  where   s   is   the  number   of    splits   required  to   yield 

a  fraction   with   no    marked  grains.       This   assumption   was    checked 

and   confirmed  by  Wright    [llj  . 

When   the   tracer   grains   had  been   counted  for   all   samples 
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in  a  given  observation,    a   cumulative   grain   count   curve  was 
plotted  and  a   rate   of    sand   drift   determined    (Appendix  II). 

As   rate   of   sand   drift   was    the   dependent   variable  in  the 
statistical  analysis   of   the  problem,    a  great  deal   of   study  and 
experimentation   was   required  for   the   establishment   of   the 
criteria  to  be  used  in   determining   this   rate.      To    eliminate 
erroneous    or   misleading  drift   rates   caused  by  anamolous   swashes, 
a         significant   percentage   of   the   total   number   of    captured 
grains   was    chosen.      In  addition,    the  region   of    maximum   slope 
on  the   cumulative   curve   indicated   that   the  greatest   volume 
of   movement   was    occurring  at   that   time.      The   criteria   developed 
°    was   Z5%  of   the   total  number   of    captured  grains   and  that   this 
percentage  lie  on  the   maximum   slope.       This    established   the  drift 
time   which  allowed  for   subsequent   calculation   of   the   drift   rate 
used  in  the  analysis.      In  addition,    the   rates    of    drift   for   five 
phi      ranges    (-1   to   V)    were   determined    (Appendix  II). 

The  final   steps   in  this   phase   were   to   determine   mean  grain 

size  and   reduce  the  data   to   a  form   compatible   with  the   data 

processing  phase    (Appendices    I,    II). 

.    Of   the  48    original   observations    the   seven   employing 

anthracene   tracers   were   discarded.      Failure   of   the  sampling 

device   resulted  in  lack   of    continuity  in  two   more   observations 

and   these   were  also   discarded.      Three  additional   observations 

were  not   used  in   the  analysis   because   their   total   grain   count 

1 
Phi    =   -log2(d),    where   d   =   grain   size   diameter    in    millimeters 
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was   less   than   10$   of   the   mean  grain   count   of   all   observations 
(  Table   I )  . 

Data   Processing 

After  the   data   were   compiled,    they   were  processed   on  a 
CDC    1604   computer.      In  order   to   provide  all  positive  inputs   to 
the   computer   the   observed   data   was    corrected   to   the  arbitrary- 
reference   levels   indicated  in   Table   I.       "Whirlpool",    a   computer 
program   for    "sorting   out"   independent   variables   by   sequential 
multiple  linear   regression,    was   utilized  for   this   procedure    [7]  . 

The  general  linear    model   may  be   expressed  by: 

Y    =   BQ+B1X1  + +BkXk+e, 

where  Y   is   the   dependent   variable,    X^ »-^k   are  t^le  Para~ 

meters    measured    (independent   variables),    Bq   is  a   fixed   constant, 

Bi, ..,Bk  are   linear   coefficients,    and   e  is   a  variable 

error   term   which  includes   data    "noise"   or   parameters   that   were 
not   considered.      The  general   mathematical   approach   will   not  be 
described  here,    since  it   is   well   presented  by   Krumbein,    Benson, 
and  Hempkins    \l~\  •      A  discussion   of   the   sum   of   squares    criterion 
will   be  included   since  it   is   a   necessary  point  in   explaining  the 
results   gleaned  from    this   procedure.      Before   proceeding,    several 
definitions   are  necessary: 

Y'   =   Y    estimate  from   regression   equations 

A 

B     =   column  vectors    of   Bs    (regression   coefficients) 

g     =   uncorrected   sums    of   cross   products    of    Y    with 


30 


vO 

o 

co 

O 

1 

O 

oJr~ 

o 

1 

CO 

m 

-4 

•O 

s 

vO 

O 

vO 

o 

CM 

CO 

O 

UA 

o 

ro 

UA 

o 

CV 

ro 

o 

CO 

UA 

rH 

CO 

UA 

ro 

vO 

o 

CO 

COUO 

CO 

•■o 

CO 

LlA 

vO 

UA 

-4 

0J 

CO 

-4 

o 

-4 

cv 

rH 

CO 

CV 

UA 

-4 

CO 

UA, 

ro 

CV 

CO 

{H 

rH 

H 

CO 

ro 

H 

o 

O 

O- 

t>- 

CO 

o 

st 

CO 

r<A 

rH 

t> 

UA 

-4 

t>- 

ro 

0J 

-4 

ro 

ro 

cv 

-4 

c- 

t> 

CO 

-4 

CV 

-4 

-4 

E-i 

r> 

sO 

c^ 

vO 

-4 

nC 

nO 

NO 

nO 

nO 

UA 

C- 

[> 

SO 

sO 

vO 

r> 

SO 

vO 

vO 

c^ 

UA 

o 

UA 

o 

uA 

UA 

no 

O 

vO 

UA 

iH 

0> 

£>- 

-4 

rH 

rH 

CV 

cv 

vO 

ua 

C- 

UA 

vO 

-4 

o 

CV 

>• 

UA 

^o 

O 

-4 

o 

MD 

UA 

co 

ro 

cr- 

!> 

ro 

O 

rH 

CO 

r^ 

ua 

CO 

NO 

ro 

cv 

cv 

t> 

-4 

CV 

cv 

O 

fr-] 

OA 

rA 

-4 

CO 

-4 

UA, 

rH 

vf) 

UA 

-41 

UA 

MD 

OA 

o 

vO 

CV 

H 

vO 

> 

rH 

rH 

UA 

-4 

O 

o 

H 

CO 

CM 

o 

co 

-4 

rH 

CO 

£>■ 

rH 

O 

CO 

vO 

c^ 

cv 

O 

CO 

ro 

CV 

rH 

CV 

O 

cv 

O 

o 

-4 

UA 

* 

O 

O 
1 

rH 
1 

-4 

rH 

rH 

cv 

cv 

rH 

CV 

rH 

rH 

o 

rH 

H 

rH 

o 

H 

rH 

1 

o 

CO 

O 

UA 

rH 

ro 

rH 

o 

1 

-4 

rH 

rH 

CV 

ITS 

3 

[*- 

o 

O 

CM 

o 

vO 

CO 

nO 

vO 

CO 

^T) 

r^ 

-; 

O 

-4 

UA 

vO 

CO 

t*- 

O 

UA 

ol 

o 

r> 

ro 

s 

UA 

UA 

CO 

o 

cv 

rH 

CO 

C  ) 

CM 

CO 

cv 

rH 

ro 

rH 

CV 

CO 

CV 

CV 

cv 

CV 

rH 

cv 

ro 

cv 

C^ 

rH 

cv 

rH 

cv 

H 

rH 

r-i 

<z 

ua 

3 

Q 

Q 

o 

O 

O 

o 

O 

O 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

O 

o 

o 

Q 

O 

CD 

O 

3 

X 

CO 

CM 

NO 

CO 

o 

-4 

rH 

rH 

o 

o 

CO 

rH 

c- 

vO 

vO 

£> 

-4 

CO 

rH 

c^ 

-4 

rH 

CV 

ro 

UA 

CO 

H 

cv 

rH 

CO 

rH 

M 

fr 

H 

rH 

1 

rH 

H 

o 

O 

o 

rH 

rH 

rH 

o 

rH 

o 

1 

o 

o 

1 

O 
1 

rH 

1 

O 

O 

1 

o 

o 

rH 

1 

O 

rH 

rH 
1 

o 

1 

O 

1 

rH 

rH 

1 

o 

O 
1 

l 

8 

Q 

o 

o 

o 

Q 

o 

Q 

o 

o 

o 

O 

R 

o 

o 

o 

O 

o 

o 

o 

o 

o 

0 

O 

o 

o 

o 

o 

o 

o 

o 

S 

0 

o 

o 

o 

f ) 

O 

o 

o 

o 

o 

C^ 

o 

C  ) 

o 

a 

( ) 

o 

o 

C) 

a 

o 

O 

C.J 

o 

o 

(' ) 

LJ 

c  > 

CD 

o 

CO 

CO 

CO 

to 

CO 

CO 

CO 

vO 

vO 

nO 

nO 

■4- 

-4 

-4 

-4- 

-4 

cv 

Oj 

cv 

cv 

o 

C) 

( 1 

o 

CO 

CO 

CO 

CO 

vO 

v£^ 

vO 

o 

Q 

rH 

rH 

rH 

r~i 

cv 

cv 

CM 

OJ 

OJ 

TA 

ro 

CO 

OA 

-4 

-J- 

-4 

-4 

-4 

-4 

-4 

-4 

u\ 

UN 

u\ 

•• 

t> 

r- 

O 

o 

CV 

-4 

c- 

o 

CO 

O 

t> 

rH 

CA 

C<A 

UA 

CV 

vO 

vO 

O 

rH 

o 

o 

O 

o 

CO 

vO 

UA 

o 

O 

O 

o 

a 

i 

CV 

cv 
I 

(N 

1 

o 

ua 

ua 

UA 

O 

o 

1 

O 
1 

UA 
1 

CO 
1 

-4 

l 

CV 

1 

1 

UA 

C^ 

nO 

i 

o 

1 

CO 

1 

o 

ro 

rH 

o 

1 

H 
1 

-4 

H 

-4 

rH 

3 

o 

CO 

1 

§ 

t> 

1 

1 

1 

E-« 

CO 

ro 

C- 

o 

co 

cv 

CO 

co 

o 

r> 

co 

■4 

o 

CO 

o 

-4 

O 

o 

o 

o 

UA 

o 

O 

CO 

o 

CV 

O 

UA 

o 

CO 

o 

> 

-4 

ua 

rH 

rH 

o 

o 

t> 

vO 

3 

o 

Pi 

rH 
CV 

O 

rH 

UA 

rH 

3 

rH 

1 

UA 

H 

PA 
rH 

o 

o 

cv 

O- 

o 
cv 

O 

rH 
CV 

■o 

CO 
rH 

UA 

vO 

o 

rH 

1 

CO 
rH 

8 

co 

ua 

O 

CM 

UA 

co 

m 

O 

UA 

ro 

o 

CO 

o 

UA 

CO 

ro 

CV 

CO 

c- 

CO 

CO 

o 

o 

UA 

UA 

MD 

UA 

vO 

o 

UA 

o 

s 

co 

CO 

ro 

UA 

-4 

CM 

ro 

-4 

no 

■4 

-4 

CV 

<r 

CO 

CO 

CO 

UA 

vO 

rH 

CO 

-4 

r-A 

UA 

CO 

ro 

UA 

-4 

UA 

UA 

nO 

UA 

g 

(V, 

£ 

rJ 

(V 

rn 

(Tri 

ni 

hJ 

(V. 

cn 

en 

tz 

H-1 

W 

Ek 

cn 

a 

tA 

K 

Cd 

J 

3 

h-1 

Crf 

Cri 

ci 

prj 

►J 

v.n 

r- 

vO 

-4 

3 

O 

UA 

O 

cv 

CO 

UA 

UA 

vO 

cv 

C^ 

rH 

CV 

CV 

CV 

U 

UN 

-4 

CV 

-4 

cv 

CO 

nJ 

H 

rH 

CV 

rH 

rH 

rH 

rH 

r-\ 

rH 

^ 

§    XI 

•H  SS 

-p 

ct3 

co 

vO 

UA 

cv 

CM 

UA 

rH 

O 

CO 

0- 

rH 

o 

UA 

CV 

CV 

UA 

ro 

UA 

UA 

rH 

o 

rH 

O 

C-- 

o 

o 

C^ 

UA 

t> 

NO 

CV 

CV 

rH 

H 

■4 

OJ 

rH 

CV 

co 

CO 

UA 

-4- 

UA 

cv 

UA 

OA 

~u- 

-4 

<t 

rH 

■4 

UA 

-4 

vO 

CO 

-4 

-4 

-4 

UA 

CV 

-4 

<f 

£    <D 

Q>  X> 

to  d 

O  23 

<< 

< 

Aj-4- 

< 

< 

■flj 

CQ 

CQ 

CQ 

CQ 

o 

o 

O 

O 

O 

Q 

Q 

C3 

a 

w 

W 

!-•-! 

w 

U^ 

fc 

fa 

fe 

CjJ 

CJ 

a 

1 

1 

cv 

1 

UA 

1 

1 

1 

CV 

1 

ro 

1 

-4 

1 

UA 

1 

rH 

CM 

CO| -4 

UA 

rH 

CV 

rO 

LT\ 

H 

CO 

-cr 

UA 

<-i 

cv 

cO-4 

cv 

-4- 

L£X 

31 


H-3 

•  UA  a" 

\  ir 

»  r-i 

c-  o- 

o 

« 

« 

i 

i 

4 

1 

1 

1 

• 

!>H 

c 

*o 

OJ 

f- 

*  [> 

r- 

-o 

vC 

o 

-3 

jy, 

CO 

-4 

<f 

CV 

ir\j^i 

o 

* 

• 

• 

• 

• 

• 

• 

• 

• 

Eh 

vO 

f- 

LO 

vO 

vO 

o 

vO 

SO 

o 

• 

-3 

m 

L^ 

to 

rH 

r-l 

cv 

C" 

o 

> 

• 

• 

« 

* 

% 

• 

• 

• 

« 

CO 

O 

o 

U"> 

01 

O 

o 

cv 

CO 

o 

£ 

fc> 

H 

rH 

rH 

£ 

0) 
O 

H 

t- 

H 

cv 

O 

C> 

-4 

c^ 

o 

• 

* 

♦ 

• 

* 

* 

* 

• 

• 

x: 

Eh 

cv 

cv 

O 

r^ 

-4 

-J 

o 

rH 

irv 

CO 

X 

1 

o 

J 

TJ 

o 

m 

C"- 

u^ 

s0 

CO 

vO 

O 

R 

S 

o 

R 

H 

R 

cv 

rH 

rH 

rH 

r-l 

< 

o 

ur\ 

• 

* 

•     • 

• 

• 

• 

• 

• 

Q 

O 

o 

o  o 

o 

o 

o 

o 

g 

o 

« 

. 

g 

R 
O 

•"—N 

'd 

E-i 

-d 

<L> 

*r 

H 

H 

CV 

o 

L/\ 

sO 

-d- 

cv 

O    r- 

o 

s 

3 

• 

• 

• 

* 

• 

• 

• 

• 

Eh 

CO 

. 

5 

rH 

r-l 

o 

o 

o 

r-l 

rH 

o 

O 

o 

1 

1 

1 

1 

Q 

PL, 

H 

o 

-p 

H 

s 

C 

CO 

CO 
rH 

o 

o 

R 

o 

nn 

J> 

»^-^ 

o 

o 

o 

o 

o 

o 

o 

o 

^3 

o< 

o 

o 

o 

o 

o 

o 

o 

o 

^r^- 

o 

JL, 

"m 

o 

NO 

<t 

-tf 

-4 

-4 

cv 

o 

CO 

g 

R 

• 

o 

Q 

LTv 

vO 

•o 

vO 

vO 

r", 

-4 

-J 

o 

«H 

9 

a 

,  i 

@ 

o 

,o 

o 

o 

a: 

CM 

CV 

~J- 

C-- 

vO 

r—t 

r-H 

E-t 

o 

E-t 

'i-s  -  . 

• 

• 

• 

• 

• 

• 

• 

• 

o 

& 

• 

H^ 

q 

o- 

H 

\0 

CV 

ts 

r-l 

o 

c^ 

fiy 

o 

3 

^ 

1 

CV 

1 

1 

1 

1 

H 

6 

o 

O 

l 

B 

o 

H 

■              r» 

li( 

pa 

+J 

c^\ 

H 

o 

o 

o 

CV 

u^ 

^0 

S 

O 

• 

• 

* 

• 

* 

* 

« 

• 

• 

s 

S 

t>- 

NO 

rH 

o 

o 

CO 

CO 

CO 

8 

o 

> 

H 

«H 

rrj 

o 

o 

CM 

co 

cv 

o 

o 

o 

O 

» 

• 

• 

• 

• 

• 

• 

• 

• 

0) 

»■— . 

ir\ 

-3 

cv 

-tf 

r> 

o- 

U"\ 

rr\ 

o 

o 

•H 

-P 

a 
o 

13 

•3 

R 

pj 

Pi 

R 

R 

>—  \ 

s 

Cti 

o 

-* 

fr- 

en 

m 

o 

r- 

cv 

irv 

«H 

ctf 

<f 

•H 

U 
QJ 

> 

CO 

U 

?H 

3 

C 

O 

O 

cv 

e<y 

Cxi 

cv 

CO 

to 

O 

o 

fl 

o 

x> 

• 

• 

• 

* 

* 

* 

* 

• 

• 

«H 

•H 

K 

O 

Lf\ 

m 

vO 

vO 

ir\ 

cv 

rr> 

o 

CO 

■P 

1C( 

?J 

ttJ 

CO 

fl 

£ 

CO 

r3 

d 

*«-^* 

T 

T 

¥ 

co   G 

o  s 

1     ' 

!  .^ 

i 
:V 

I 

1 

•71 

i 

% 

* 

32 


successive  X's 
SSy  =  sum   of   squares    of    Y 
The  first   step   of   the   procedure  is   to   compute   SSy   ^y 
the   equation: 

SSY    =    £(Y-Y)2  (2) 

by  using  the   computing  formulas 

ssY  =  £y2-y]Ty  ,  (3) 

Since   the   mean  value   of   the   observed   Y   values,    Y,    is   the   same 
as   the   mean  Y1   of   the   computed  Y'   values,    then  SSyi    can  be 
computed  using: 

ssYi  =]Ty,2-y^y  (4) 

where 

ZY'2   =  i.T£  (5) 

If  Yy'2   =  ^~Y2,    this    would  indicate   that   the  variability   of   Y   and 
Y'   were   the   same   from    (3)    and    (4)» 

The   criterion    for      ranking   the   parameters   in   their   order 
of   importance  is   based   on  the   degree   to    which  any   single   para- 
meter  or   combination   of   parameters    "accounts   for"   the 
variability   of   Y.      In  other   words,    if   Y'   for    each   set   of   para- 
meters  was   exactly   equal   to    each   Y    observed,    then   100$   of   the 
variability   of   Y    would  be    "accounted  for".      The  percentage   of 
the  variability   each  parameter    or   combination  of   parameters 
"accounts   for"   Y   is    computed  from: 

Percent    Red*-    100(SSy,)/(SSY)  (6) 

This    quantity         is        informally    called   the    "percentage   SS   reduction 
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of   Y"   and  the   machine   output   lists   this    quantity  for   all   stages 
of   the  analysis.      The   reduction  in   the   SSy   is   a   measure   of   the 
mathematical  association  between  variables    (here  the  linear 
association)    and  it  is   not   necessarily  a   measure   of   the  physical 
relationship. 

To   evaluate   the  relative  importance   of   any  parameter   or 
combination  of   parameters    on   this   least-square  basis,    the 
scanning   routine  involves    the   selection   of   the   strongest   single 
parameter,    the   strongest   pair   of   parameters,    and   continues 
until   the  strongest   combinations   are   examined. 

To   illustrate  the  foregoing   statement,    a  fictitious    machine 
output   for   three  parameters   is   given  in   Table  2. 

TABLE  2t       FICTIOUS    EXAMPLE     OF    THREE   X'S 

X'S  NUMBER    OF   X'S 

ONE   AT   A   TIME 
1 
2 

3 

TWO  AT  A   TIME 
1,2 
1,3 
2,3 

THREE  AT   A   TIME 
1,2,3  75.0 

The   strongest   parameter   is    X->    with   a    "percent   SS   reduction" 

of   50$,    the  strongest   pair    of   parameters   is   X.,    and  X^   with   a 

"percent   SS   reduction"    of   70$.      From   this,    the  net   conbribution 

of   X-l,    in   combination   with  X?,    is    70$-50$   =   20$.      This    difference 

from   the  value   of   X^   alone,    25$,    essentially  indicates   some 
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PERCENT 

SS 

REDUCTION. 

25- 

,0 

50, 

,0 

12, 

-5 

70, 

,0 

50, 

,0 

25, 

,0 

redundancy   exists   between  X^   and  Xg. 

In  this   example  Xo    is   obviously  weakest,    and  its   contribution 
In  the  presence  of  both  X^   and  X^  is   75^-70$  =   5%*      Comparison 
with  X~    taken  alone  also   shows   that   some   data   redundancy  is 
present. 

In  this    example,    the   difference  between  the  parameters 
is      sufficiently   large   so   that   ranking   them   as   to   importance  is 
relatively   easy,    but  in   other   real  problems,    the  differences    may 
be   small   so   that   ranking  is   at   least  partially   subjective. 

If   non-linear   relations   are   suspected  or  known,    conventional 
transformations   of    the  parameters    may  be   made  before  using 
the  program.      Transformations    such  as   log  X,    X2,    e~x,  etc., 
can  be  used  to  help    "reduce  the  SS".      A  logarithmic   trans- 
formation of  distance  from   the   canyon  was   used  in  this   study. 

Several  items   should  be  noted  here  before  proceeding   to 
the  results  and  conclusions,    they  are: 

(1)  A  small    "percentage   SS   reduction"   indicates    either 
erroneous   data   or   omission   of   important  variables; 

(2)  Several   combinations    of   parameters    can  be   made  which 
will   yield  different    "percent    SS   reduction"   for   the   same  variable; 
and, 

(3)  Ranking  of   the  parameters   can  be   different   depending 
on  subjectiveness    of   the  investigator. 

Items    (2)    and    (3)    result   from    the  fact   that   the   strongest 
combination  which  fits   with   the   previous    combination   might   not 
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be  the  first   strongest,    but    may  be  the   second  or    third,    whereas 
the  first   strongest   will  fit  in  a   different   series.      All   variables 
eventually   make  an  appearance,    although  their   orders   in  various 
series   are  different,    with   correspondingly   different  values.      A 
series   is   like  that   described   in  the  fictitious    example  and   Table    2, 
iie.',.       a   series        for   this    example  is   X-^;    XifX^;    XifX^Xo. 

The   distinct   advantage   of   this   procedure  over   other   step- 
wise multiple  linear   regression  programs,    i.e.,    the  BMD    series 
developed  by   UCLA  L^J ,    is   the   capability  of   looking  at   the 
many   combinations   of   parameters.      Other   programs   create  a 
new    combination,    e.g.,    four  at  a   time  from   three  at  a   time, 
by  adding  a   new   parameter   that   gives   the   "best   goodness   of 
fit."      As   can  be   seen  from   the   comparison  between   "Whirlpool" 
(Table  3)    and  a   typical  BMD   output    (Table  4)    the  ranking   may 
be  different   and   erroneous   for   the  BMD   approach. 
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TABLE  4;   A  TYPICAL  BMP  TYPE  OUTPUT  FOR  A 
STEPWISE  MULTIPLE  LINEAR  REGRESSION  PROGRAM. 


RELATIVE  RANKING 
WITH  OTHER 
COMBINATIONS   X'S 


PERCENT  SS 
REDUCTION 


1 
1 
2 
4 
3 
6 
1 
1 


D, 


DC. 

>Hb 

•Dd 

rHb, 

DC! 

rHb, 

DCi 

rHb, 

DCi 

»Hb, 

DC. 

-Hb, 

Poi 

•Hb, 

D50 

■D50»a»HT 

D50»a»HT»VLONG 

D50»a»HT'VLONG,T 
D50,a»HT»VLONG»T*VON 


15.98 
27.10 

31.33 

36.14 

41.25 
43.36 
46.67 
48.43 
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5.       RESULTS 

The  results    of    the  investigation  are  presented  in  two 
parts   dealing  with  the   subproblems. 

Direction   of    Sand  Movement 
Table   1   is  a   tabulation  of   the   data  gathered  during   the 
various   observation  periods.      A  perusal   of   this   table  reveals 
that   all  36   observations   used  for   final  analysis   indicated  positive 
sand  movement  toward  the  canyon.      In  no   case  was   there  move- 
ment  in  the   opposite   direction. 

Rate   of   Sand  Movement 

The  results   for   this   particular   problem   were  the  outputs 
of   the    "Whirlpool"   program.       Table   5   presents   the  five  best 
"percent   SS   reduction"   for    eight   combinations    of    eleven 
independent  parameters.      The   maximum   total   "percent   SS 
reduction"  for  this   set   of   data  was   51.54$    (linear)    and  58.02$ 
(non-linear)    for   all   eleven  parameters    considered  together.      As 
can  be   seen,  from   Table  3,  the   maximum    "percent   SS   reduction" 
for   eight   linear  independent   variables   was   48 -43$   which  was 
considered  sufficient   since  this  accounted  for  93-97$  of   the 
original   "percent  SS  reduction".      Table  3A  shows   that   the 
maximum    "percent  SS   reduction"   for  the  non-linear   case  with 
six  independent  variables    was  54. 67$  which   accounted  for  94-21$ 
of   the  original    "percent   SS   reduction". 

Table  3   lists   the   contributions   by   each   parameter   among 
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TABLE  5 :  THE  STRONGEST  PERCENT  REDUCTION  IN  SAND  DRIFT  SUM  OF  SQUARES 
ATTRIBUTABLE  TO  EACH  OF  SEVERAL  COMBINATIONS  OF  ELEVEN 
INDEPENDENT  VARIABLES  (TOTAL  PERCENT  REDUCTION,  ALL  X»S  51.54) 

X»S  NUMBER  OF  X»S  PERCENT  SS 
'_                                                                                REDUCTION 


ONE  AT  A  TIME 


DC 

CL 

Hx 

VON 

M 

TWO  AT  A  TIME 

Dc,Hb 

DC>* 

Dc,Ht 

Dg,M 

DC»vON 

THREE  AT  A  TIM 

Dc,HT,a 

Dc»Hb>D50 

Dc»Hb>a 

Dc>Hb,V0N 

DC,Hb,T 

FOUR  AT  A  TIME 

DC,HT,a,T 

Dc#HT,a,D5Q 

. 

Dc,Hb,V0N,T 

Dc*Hb>a,D5Q 

DC*Hb»HT'D50 

FIVE  AT  A  TIME 

Dc,HT,a,T,VLONG 

Dc*HT,a,T,V0N 

Dc,HT,a,D^0,Hb 

Dc,a,T,Hb,V0N 

DCa>T»HT'D50 

SIX  AT  A  TIME 

Dc#a,T,HT,Hb,VL0NG 

Dc,a,T,HT,V0N,VL0NG 

Dc»a»T*HT»D50#vLONG 

DC'a'T'HT'CL,VLONG 

^C'a'   '^T'^b'   ON 

15.98 

10.51 

3.52 

2.60 

1.63 

27.10 
25.61 
25.27 
20.47 
19.57 

32.46 
31.33 
31.09 
31.05 
29.95 

37.93 
37.12 
36.75 
36.14 
35.69 

42.65 
41.28 
41.25 

40.76 
40.21 

44.32 
44.21 
43.60 
43.53 
43.38 
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TABLE  5  (continued) 
X«S 


NUMBER  OF  X*S 


PERCENT  SS 
REDUCTION 


Dc,a,T,HT,Hb,D50,VLOrjG 

DC*a#T*HT»Hb»vON»VLONG 
Dc,a,T,HT,Hb,VON,D5o 

DC'a'T»HT»CL»D50»VLONG 
Dc,a,T,HT,CL>Vs,VQN 

DC>a#T*HT#Hb#D50»VON*VLONG 
DQjajTjHrpjH^DjQfVQj^Vg 

Dc,a,T,UT,Hb,V0N,VL0NG,Vs 
Dc,a,T,HT,CL,V0N#VL0NG,Vs 

Dc#a#T>HT#Hb*D50'VLONG'VS 


SEVEN  AT  A  TIME 


EIGHT  AT  A  TIME 


46.67 
46.52 

46.35 
44.97 
44.94 

48.43 
48.34 
47.70 

47.43 

47.33 


43 


the   eight   strongest   independent   variables   taken  jointly.      It 
should  be  noted  that   levels    two,    six,    and   seven,    are  the   second 
strongest   combinations,    since   the  first   strongest   did  not   mesh 
into   the   scheme.      It   is    evident   that   distance  from   the   canyon 
was   the   strongest   contributor   with  a    "percent   SS   reduction" 
of   15.98$.      Even  though  longshore   current  and   mean  slope 
appeared  in  the  five   strongest  individual   contributors,    they   did 
not  appear  in  the  final   analysis.       Table   3-A.  lists   the   contributions 
by   each   parameter   among   the   six  strongest'  independent   variables 
taken  jointly.      As    contrasted   to   the   linear    case,    all   levels   were 
the   strongest   combinations.       The   strongest   contributor   was 
again   distance  from    the   canyon   with   a    "percent   SS    reduction" 
of    17.86$. 

Table  6   shows    the  range  and   mean  value   of    sand   drift   at 
each   station  over   the   period   of   investigation.       The   trend    of 
the  mean  rates   was   inversely   proportional   to   the   distance  from 
the   canyon.      The  regression   equation  that   resulted  from   the 
linear   analysis: 

Y   =   14.403+0.i+96Hb+0.109a+0.147M-0.069Von-0.069Vlo 

-0 .  001Dc+0  .  514C1-8 .364Dc0+0  .  524Ht+0 .  172Vs-0 .921T    (7) 
yielded   the  following  answers    when  the   data   from   2+D,    2E,    and 
5F    (corrected   to   reference  levels,    Table   2)    were  inserted: 


kk 


Sample           Observed            Computed  Percent 

Number Drift Drift Residuals        Accuracy 


4D 

11.1 

9.9 

2E 

6.7 

5.1 

5F 

6.7 

5.7 

1.2 

89.2 

1.6 

76.1 

1.0 

85.1 

The  average  accuracy  of   these  tests   on  the  regression  coeffic- 
ients  was   83. 5% •      It  was   previously  noted  that  these  three 
observations   were   excluded  from   the  data   processing  due   to 
insufficient  grain   count. 
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TABLE  6 :  RANGE  AND  MEAN  RATE 
OF  SAND  DRIFT  :AT  EACH  STATION 


STATION  • 

DISTANCE  FROM 
CANYON 

RANGE  OF 
SAND  DRIFT 

MEAN 
RATE 

A 

800 

3. 3-16 ,7  FT/MEN 

6.4FT/MIN 

B 

1600 

5.0-8.4  FT/MIN 

7.0FT/MIN 

C 

2400 

2.6-6.7  FT/MIN 

4.9FT/MIN 

D 

3200 

0.8-4.2  FT/MIN 

3.2FT/MIN 

E 

4000 

1.9-7.1  FT/MIN 

3.5FT/MIN 

F 

4800 

5.0-10.0FT/MIN 

7.4FT/MIN 

G 

5600 

2.5-3.3  FT/MIN 

3.1FT/MIN 

H 

6400 

0.0-7.5  FT/MIN 

3. 3  FT/MIN 

6.       CONCLUSIONS 

The   sampling   technique   was   designed  to   provide  distinction 
between   dispersion  and   drift   in   either   direction.      Sands   injected 
toward  the  canyon    (orange)    were  never   found,    while   those 
injected  away    (green)    were  found  in  all   observations.      These 
results   indicate  that   dispersion  and  drift   away  from   the   canyon 
did  not   occur   during  this   study   period.      It   can  be   concluded 
then,    that   during   the  winter/early-spring   season,    within  the 
area  investigated,    the   swash-zone-  sand   movement   is   directed 
toward   the   canyon. 

The  parameters   chosen  have  a   definite,    though  not 
completely  understood,    effect   on  sand  drift.      The  assumptions 
made  were   carefully   considered  and  are  believed   to  be   correct 
as  far  as  beach  sand  and  ocean  interactions  are  known.      With- 
out  more  concrete  knowledge   of   the  interactions,    the  only 
conclusions   which  can  be  drawn  are  statistical  in  nature. 
Consequently,    it  is   felt  that  a    "Whirlpool"   type  program   is   a 
superb    method  for   ranking   the  important  parameters. 

Since  in  the  linear  and  non-linear   cases   there  was  a 
maximum   of   58.02$    "percent   SS    reduction",    this   is   an   excellent 
indication  that   imprecise   techniques   in  the   measurements   led   to 
data    "noise".      The   methods   described  in  Appendix    I    involved 
many  approximations.       Exact    measurements   of   these  para- 
meters,   however,    would  have   required   more   elaborate  instru- 
ments  than  were  available. 
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Another   source  of   data    "noise"    could  have  resulted  from 
a  poor   choice   of   criteria  in  determining   rate   of  drift.       Even 
though  the   "25^/steep   slope"   criterion  gave  the  best  results   of 
those  evaluated,    there  may  have  been  a   more  optimum  selection. 

Of   the  parameters   investigated,    there  are  two   that    may 
not  be  indicative   of   the   energy   relationships   which  they  were 
chosen  to   represent.      There  are  two    observable   components 
of   swash-velocity,    onshore  and  longshore.      Only  the   onshore   one 
was  used  in  this   study.      A  more  desirable  parameter    may  have 
been  the  longshore  component,    but  this   value  was  not  used 
because  no   observation  of   the  nearshore  angle  of  approach  was 
made.      The  value  used  for   mean  grain  size  was  a  function  of 
the  drift   rate  time.       Since   the  latter   was   based   on  arbitrary 
criteria,    the   mean  grain   size   would   similarily  contribute  to   data 
"noise".      Therefore,    the   longshore   component    of   swash-velocity 
might   replace   onshore   swash-velocity  and  .mean   grain   size;  as  'being 
more  representative   of   the   energy   distribution  on  and  along 

the   swash-zone. 

The  final   parameter   suspected   of   generating   data    "noise" 
was   longshore   currents.      These  values   were  the  results    of 
short-term   observations    (two    minutes    or  less).      More  represent- 
ative  values    might   have  been   obtained  by   measuring   long-term 
mean   currents. 

Another   factor   contributing  to   this   relatively   small    "percent 
reduction"   was   the  limited  interval   of   time   over   which   the   study 
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was   conducted.      The   study   did  not  have  a  sufficient  time   span 
for  the  variability  of   the  parameters   to   reach   their   maximum 
possible  range  of   values. 

The  large   value   of    "percent   SS   reduction"   for   distance 
from   the   canyon   may  be  partially  attributed  to   the  fact   that 
there   were   only   eight  variations   of   the   distance  available, 
causing  a  possible  biasing   effect.      More   random   values   of  the 
distance  input   may  have  yielded  less   promising   results. 

Ignoring   the  limitations    of   short-term   observations   of   long- 
shore  current,    the   low   individual    "net   percent   SS   reduction" 
obtained  for   this   parameter   can  be  interpreted  to   mean  that 
longshore  current  does   not   significantly   contribute   to   the   swash- 
zone  sand   movement. 

In  view   of  the  preceding  facts,    no   concrete  conclusions 
can  be  drawn  concerning   the  relative  importance  of  the  observed 
parameters   in  affecting  the  rate   of   drift.      In  light   of   earlier 
energy   considerations,    it   is   not   unreasonable,    however,    to 
assume  that   the  relative  ranking   of   the   parameters   in  the  non- 
linear  case  is   indeed  indicative   of   the   real   situation.      Possible 
exceptions  to  this  ranking  include  the  apparently  negligible 
contribution  by   mean  slope  and   longshore   current. 

This   study   was   not   completely  restricted  to    evaluating 
only  the   statistical   results.      Five  individual  phi   range  rates 
were   determined  but  not   used  in   the  statistical  analysis. 
Inspection  of   these   rates    revealed   that   with   increasing  grain 
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size,    the  rates   decreased   proportionately   to   the  point   where 
no   movement   existed  in   the  largest  phi   range    (Table  7).      This 
fact   implies   that   larger   grains    either    exhibit  no   longshore   move- 
ment  or   that   the  rates   were   too   small   to  be   detected   during 
the   observational  period. 

Another   interesting  aspect   of   this   study   was   the  highly 
satisfactory   performance  of   the   sampling  device.      It   is   an 
improvement   over   other   sampling  techniques   in   that   a   single 
observation  ■will  provide   conclusive   evidence   of   direction  of    sand 
movement   or   dispersion. 

Finally,    angle   of  approach,    wind  alongshore,    and  longshore 
current   should,    by   classical   concepts,    govern  the  direction  and 
somewhat   influence  the  rate   of   sand   movement,    providing  the 
study  area  has  a  straight  beach  with  relatively  parallel  offshore 
contours    [_10 J  .      It   is   obvious   from     the   data   in  Table   1   that 
in   many  of   the   observed   cases   these  relationships   did  not   hold 
for   the   swash-zone.       Therefore,    it   can  be   concluded  that   either 
the   swash-zone   movement   is   affected  differently   than  the    net 
littoral   movement,    or   that   the   canyon  has   an   overriding   effect 
on  the   movement   of   sand  in  the  swash-.zone. 
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TABLE  7 :  SAND  DRIFT  RATES  FOR  EACH 
OBSERVATION  FOR  PHI  RANGES  -1  THROUGH  4 


OBSERVATION 

PHI  RANGES 

NUMBER 

-1  TO  0 

0  TO  1 

1  TO  2 

2  TO  3 

3  TO  4 

1-A 

1.4 

2.5 

3.6 

3.8 

3.1 

2-A 

0.0 

2.0 

2.0 

'   5.6 

3.3 

3-A 

0.0 

3.2 

3.8 

3.3. 

3.3 

4-A 

2.5 

2.9 

6.7 

3.3 

5.0 

5-A 

0.0 

0.0 

4.5 

5.0 

16.7 

6-A 

0.0 

2.0 

3.8 

2.0 

3.6 

7-A 

0.0 

2.3 

2.4 

2.5 

2.5 

2-B 

1.8 

4.2 

5.0 

4.5 

4.2 

3-B 

0.0 

1.9 

2.0 

7.1 

2.0 

4-B 

2.8 

6.3 

5.6 

5.6 

5-B 

0.0 

3.3 

3.3 

16.7 

25.0 

1-C 

2.5 

1.7 

3.8 

4.5 

8.3 

2-C 

0.0 

6.3 

6.7 

6.3 

6.7 

3-C 

0.0 

5.0 

5.0 

5.6 

6.3 

4-C 

0.0 

0.0 

3.3 

3.6- 

8.3 

5-C 

Q.O 

2.3 

2.5 

2.5 

12.5 

1-D 

0.0 

1.4 

2.5 

4.2 

5.0 

2-D 

0.0 

2.8 

3.8 

4.5 

12.5 

3-D 

0.0 

0.0 

5.0 

7.1 

10.0 

5-D 

1.7 

2.0 

3.6 

3.3 

16.7 

1-E 

0.0 

2.6 

3.3. 

8.3 

12.5 

3-E 

0.0 

1.8 

1.9 

1.8 

2.0 

4-E 

0.0 

3.1 

3.3 

8.3 

12.5 

5-E 

0.0 

1.8 

2.5 

3.3 
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7.   FUTURE  STUDY 

Based   on  the  assumptions    made  in   this   thesis,    extending 
the   study  to   include   sampling:    at   greater    distances    on  both 
sides    of   the  canyon;    in  the   canyon;    in  the  surf    zone;    and 
through  a  vertical   profile,    would  be   highly  informative. 
Additionally,    the  study  period   should  be   lengthened,    and   more 
precise   techniques    should  be   used   to   measure    the   desired 
parameters. 


53 


8.      ACKNOWLEDGMENTS 

We   extend   our  appreciation   to   our  advisor,    Assistant 

Professor   Victor   T.    Neal,    Department   of   Meterology  and 

Oceanography,    United   States    Naval   Postgraduate   School,    for 

his   assistance   throughout   the   period   of   this    study.       The  authors 

are  also   indebted  to   the   United   States  Army   Infantry   Training 

Command,    Fort   Ord,    California  for   the   use   of   a   «Jeep   vehicle 

during   the  field   phase   of   this    study,    and  in  particular,    the 

personnel   of   the  Weapons   and   Equipment  Motor    Pool   for   their 

excellent   cooperation  and  assistance  in  providing   the  vehicle  for 

daily  use.      Additionally,    thanks   are   due  to   Mrs.    Davis   and 

Mrs.    Harper, who    donated   their   time,    ovens,    and   patience   while 

drying   the  many   sand   samples. 


55 


BIBLIOGRAPHY 

1.  Chamberlain,      T.K.,    "Mass   Transport   of    Sediments   in  the 

Heads   of   Scripps    Submarine   Canyon   California", 
Papers   in  Marine  Geology,    pp.    42-6V,    Macmillan 
Company,    New    York,    1964. 

2.  Dixon,    W.vJ.,    ed.  ,    BMP,    Biomedical   Computer    Programs,- 

Health   Sciences    Computing   Facility,    Department   of 
Preventative  Medicine  and   Public   Health,    School   of 
Medicine,    University   of    California,    Los   Angeles,    196i+. 

3.  Galvin,  ,C  .  «J .  ,    A   Selected  Bibliography  and   Brief   Review 

of   the   Theory  and    Use   of   Tracers    in   Sediment 
Transport    Studies,      CERC   Bulletin   Summary   Report 
Vol.    I,    U.S.    Army   Coastal   Engineering   Research 
Center,    Washington,    D.C.,    196i+,    pp.    5-12. 

4.  Harrison,    W.    and   Krumbein,    W.C.,    Interaction   of   the 

Beach-Ocean  Atmosphere   System   at   Virginia  Beach, 
Virginia,    U.S.    Army   Coastal    Engineering   Research 
Center,    Tech.    Memo.    No.    7,    Washington,    D.C.,    I964 

5.  Hohenstein,    C.G.,    Jaeger,    JT.W.  ,    and   Jones,    D.L., 

A   Study   of   Marked   Sand  Movement   on   Del  Monte 
Beachj    Monterey  Bay,    California,    M.S.    Thesis, 
U.S.    Naval   Postgraduate   School,    Monterey,    Calif- 
ornia,   I965. 

6.  Krumbein,    W.C.,    The  Analysis    of   Observational    Data 

from   Natural   Beaches,    Dept.    of  Army,    Corps    of 
Engineers   Beach    Erosion  Board,    Tech.    Memo.    No. 
130,    1961. 

7.  Krumbein,    W.C.,    Benson,    B.T.,    and   Hempkins,    W.B., 

Whirlpool,    A   Computer    Program    for    "Sorting   Out" 
Independent  Variables   by   Sequential  Multiple  Linear 
Regression,    Tech.    Report    No.    li+,    ONR    Task    No. 
389-135 »    Office   of   Naval   Research,    Geography 
Branch,    Northwestern   University,    I964. 

8.  Sheppard,    F.P.,    Mass   Movement   in   Submarine   Canyon 

Heads,    Trans,    Amer.    Geophys.    Union,    32,    3»    PP» 
405-418,    June   1951. 

9.  Sheppard,    P.P.,    and   Inman,    D.L.,    Nearshore  Water 

Circulation   Related   to   Bottom    Topography  and  Wave 
Refraction,    Trans.    Amer.    Geophys.    Union,    31»    3» 
pp.    196-212.    April    1950. 

57 


10.  Wiegel,    R.L.,    Oceanographic    Engineering,    New   Jerseys 

Prentice-Hall   Inc.,    19 64. 

11.  Wright,    F.F.,    The    Development   and  Application   of   a 

Fluorescent   Marking   Technique  for    Tracing   Sand 
Movements   on  Beaches,    U.S.    Navy,    ONR,    Geography 
Branch,    Project   N.R.    388-057,    Tech.    Rept.    2, 
(Columbia   University,    Dept.  of   Geology),    1962. 

12.  Yasso,    W.E.,    Use   of   Fluorescent   Tracers    to   Determine 

Foreshore    Sediment    Transport,    Sandy    Hook,    New 
Jersey,     U.S.    Navy   ONR,    Geography  Branch,    Project 
N.R.    388-057,    Contract   Nonr.    266(68),    Tech„    Rept. 
6,    (Columbia   University,    Dept.    of   Geology),    196i+» 


58 


APPENDIX  I 

DESCRIPTION  AND  MEASUREMENT 
OF  THE  STUDIED  VARIABLES 


Hb 

» 

Height    of   the  Breaking  Waves 

A   man   was   placed  at   the   mean   water   level   with   a   stadia 
rod   while   another    sighted   across    the    mean  breaker   level, 
recording   the  height   on  the   stadia   rod.       The   average   of   these 
measurements    taken   during   the   observation  period   was   used 
for   the   value   of   breaker   height. 


Angle   of   Approach 

A  Brunton   Compass    was   utilized  in   measuring   the  angle 
of   approach.      One    man   pointed   the   compass   normal   to   the 
approaching   wave   crests   while   another   man  lined   up  the 
beach  face.      At   least   two    measurements    were   taken  during 
each   observational   period.       The  final  value '  was   taken  to  be 
the  average   of   these   measurements. 

M 

Mean   Slope   of    the   Swash-zone 

A  Brunton   Compass    was   used   to    measure   the  beach   slope 
in  the   mid   swash-zone. 

VON    and   VLONG 
Onshore   and  Longshore  Wind 

A  standard   Navy   wind   measuring   set,    AN/PMQ-3C,    was 
used   to    measure   wind   speed  and   direction.      Three   values    were 
taken  per   observational   period.       These  values,  were  averaged 
and  broken  down  into   onshore  and   longshore   components. 
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Dc 

Distance  from   the   Canyon 

Distances   from   the   canyon  -were   surveyed   with   a   transit 
and   stadia  rod  using   the  Moss   Landing   Pier  as   a  base  point. 

CL 
Longshore   Current 

Three  to  five  drift  bottle  measurements  were  made  over 
a  100  foot  run  during  each  observational  period.  These  values 
were  averaged  and   the  result   designated  as   longshore   current. 

D50 

Mean   Grain   Size 

Mean  grain   size  was  d'eter  mined  from   a   cumulative  percent 
curve  of    the  fluorescent   tracer   grains    captured  up   to   and 
including  the  time   designated   as    the   Drift   Rate   Time. 

ilrp 

Height   of   the   Tide 

Height   of    the   tide   was    calculated  from    standard   tide 
tables   for   the    time  halfway   through   each   observational   period. 


Period   of   the  Waves 

Wave  period   was    measured  at   the   start   of,    halfway 
through,    and  at   the   end   of    each   sampling  day.      The   periods 
of   30    waves   passing   the   end   of   the  Moss  Landing   Pier   were 
averaged  for    each   wave   measurement.      The   resulting  values 
were  applied   to   the  nearest    observational  period. 

VS 
Swash- Velocity 

Swash-velocity   was    determined  from   two    measured   values, 
swash-length  and   swash-period.       The   swash-length  was    the 
average  length  of    15   or    more  swashes.      The   swash-period   was 
the  average   of   a   series   of   full   swashes.      A  full   swash   was 
considered  to   be  one   that   was   not  noticeably  interfered   with 
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by  backwash   and   uprush.  The  average   swash-length   was 

doubled  and   then   divided  by    the   average   swash-period  to    give 
a   value  for   average   swash-vfjlocity  for   the   observational 
period. 


Rate   of    Sand   Drift 

Rate  of   sand   drift   is   discussed  in   the   Data  Analysis 
section. 
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APPENDIX    II 
Sample   Data  Analysis    Sheets 
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SAMPLE 

DATA    OBSERVATION  SHEET 

Station        D       Observation   Number        2  Date      27    Feb.     Time     1130 

Breaker   Height    (Hb)  4.  5    (ft) 

Angle  of  Approach    (a)  1R  (  °) 

Mean  Slope  of   Swash-Zone    (  M)  6.3    (  °) 

Wind   Speed  10.0    (kts) 

Wind  Onshore(V        )  9.7    (kts) 

Wind  Dongshore(  v"long)  2.6    (kts) 

Wind   Direction        345    (  ° )    Relative  to   onshore  normal 

to   beach 

Distance  from    Canyon    (Dc)  3200    (ft) 

Longshore   Current    (CL)  0.8    ( ft/sec ] 

Height   of   Tide    (HTj  1.3    (ft) 


S-wash   Period  8.9    (sees) 


•wash  Velocity    (Vs)  4.5    (ft/sec) 


Swash  Length  20    (ft) 

Wave   Period    (T)  6.3    (sees) 

Drift   Rate  Time  12.0    (min) 

Mean  Grain  Size    (Dr0)  0.25    (mm) 

Drift   Rate    (Y)  4.2    (ft/mln) 

Drift   Direction  Toward 


Remarks:      An  identical   sheet   was   prepared  for    each  observation 
made.       Direct   measurements    were  recorded  in   small   notebooks 
and   a   mean  value   obtained   prior   to   tabulation.      Subjective   deter- 
mination of   mean  grain   size,    rate,    and   direction  of    drift   concluded 
the  data   compilation  of   each   observation. 
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SAMPLE  GRAIN  COUNT  SHEET 


Sample  Number  2-D-3 


Time  After  Seeding   15 


Median 
kSize(mm) 

Phi 
Range 

No.   of 
Splits 

No.    of 
Grains 

Cumulative 
Grains 

Percent 
Total 

Cumulative 
Percent 

1.5 

-1-0 

0 

0 

0 

0 

0 

.75 

0-1 

0 

0 

0 

0 

0 

.375 

1-2 

8 

256 

328 

48,  S 

48,8 

.188 

2-3 

7 

128 

224 

33*4 

82.-2 

.094 

3-4 

5 

32 

120 

17.-8 

100.0 

iH 

Total 

416 

676 

s^SSisi 

nm 

Remarks:   A  similar  sheet  was  prepared  for  each  sample, and 
were  combined  to  plot  a  cumulative  grain  count  ourve  for 
each  observation.   The  percent  values  were  only  computed 
for  the  sample  following  the  drift  rate  time  and  were  later 
used  in  the  determination  of  mean  grain  size. 


65 


SAMPLE  CUMULATIVE  GRAIN  COUim'T  CURVE 


Sample  Number" 2-D- 


TOTAL 
GRAINS  ,00° 


I    25    30    35    40  45     50    55    60 
MINUTES 


'I 
12 


Remarks;   Similar  curves  were  prepared  for  each  Phi  range 
investigated.   Grain  counts  for  each  sample  were  plotted  vs. 
time,  and  the  rate  of  drift  was  determined  at  the  time  when:1 
25$  of  the  total  grains  lay  on  the  maximum  slope  of  the  curve, 
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SAMPLE  MEAN  GRAIN  SIZE  CURVE 


Sample  Number  2-D 


100 


.75 


%  60 


25 


PHI 


1                                              2 

\                                           t 

i/'W.2% 

He.6% 

1 

5                                       .2 

5                                         J 

25                                     .06 

25 

MM. 


'50 


Remarks:   Mean  grain  size  for  each  observation  was  chosen 
as  the  size  corresponding  to  50$  of  the  cumulative  per- 
centage of  captured. tracers. 
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text.   The  assignment  of  links,  rnles,  snd  weights  is  optional. 
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